Nonlinear stochastic biasing of halos: 
Analysis of cosmological A^-body simulations and perturbation theories 
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It is crucial to understand and model a behavior of galaxy biasing for future ambitious galaxy redshift sur- 
veys. Using 40 large cosmological A'-body simulations for a standard ACDM cosmology, we study the cross- 
con'elation coefficient between matter and halo density field, which is an indicator of stochasticity of bias, over a 
wide redshift range < z < 3. The cross-correlation coefficient is important to extract information on the matter 
density field, e.g., by combining galaxy clustering and galaxy-galaxy-lensing measurements. We compare the 
simulation results with integrated perturbation theory (iPT) proposed by one of the present authors and standard 
perturbation theory (SPT) combined with a phenomenological model of local bias. The cross-correlation coef- 
ficient derived from the iPT agrees with A'-body simulation results down to r ~ 15 (10) /7"'Mpc within 0.5 (1.0) 
% for all redshifts and halo masses we consider. The SPT with local bias does not explain complicated behav- 
iors on quasi-linear scales at low redshifts, while roughly reproduces general behavior of the cross-correlation 
coefficient on fully nonlinear scales. The iPT is powerful to predict the cross-correlation coefficient down to 
quasi-linear regimes with a high precision. 
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I. INTRODUCTION 



In the standard cosmological model, known as the ACDM 
model, the energy density is dominated by mysterious compo- 
nents called as dark matter and dark energy. The correlation 
function of dark matter and its Fourier counterpart, the power 
spectrum, contain a wealth of information that can be used 
to determine, e.g., the dark matter, dark energy, and neutrino 
masses. Thus, it is very important to exploit these quantities 
in the large-scale structure of the universe, which is a pillar 
of modern observational cosmology. However, it needs to 
be investigated how to take account of galaxy biasing. Ob- 
servable galaxies are biased relative to the underlying matter 
density field. The galaxy biasing is affected by nonlinear ef- 
fects and is scale dependent in general. Such nonlinear effects 
inipose a serious problem in analyzing galaxy surveys [e.g., 
[ll-H]. Upcoming galaxy surveys such as BigBOSS' Jsll, Eu- 
clid la], Subaru Prime Focus Spectrograph (PFS)^ ||7|l and the 
Wide-Field Infrared Survey Telescope (WFIRST)^ require the 
understanding of galaxy biasing with high precision and thus 
theoretically precise description of the galaxy biasing is a cru- 
cial issue. 

Most of the direct studies of the clustering of matter on 
cosmological scales rely on shear-shear weak lensing, but it 
is also possible to extract information on the matter cluster- 
ing by combining galaxy clustering and galaxy-galaxy lensing 
measurements [e.g., S]. To achieve this, one has to precisely 
know the relation between the distribution of galaxies and the 
distribution of matter An important property of the relation 
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is often characterized by a cross-correlation coefficient. The 
cross-correlation coefficient is a characteristic parameter of 
stochasticity ||9]. Since galaxies are expected to form in dark 
matter halos in modern models of galaxy formation, under- 
standing and modeling the clustering properties of the halos 
play an important role and are crucial first step in modeling 
galaxy biasing. 

In this work, we examine how well known models of halo 
clustering reproduce the cross-correlation coefficient between 
matter and halo density fields obtained from A^-body sim- 
ulations. We consider two models of nonlinear bias: the 
inte grat ed Perturbation Theory (iPT) developed by Matsub- 
ara lUOl] which naturally incorporates the halo bias, redshift- 
space distortions, nonlocal Lagrangian bias and primordial 
non-Gaussianity in a formalism of perturbation theory with a 
resummation technique based on the Lagrangian picture (see 
also lUU \VM ). and the standard perturbation theory (SPT) 
combined with phenomenological model of local bias, which 
leads to non-trivial renormaUzations of the leading order bias 
parameter 11 1311 . A significant advantage of the iPT is that it is 
simpler and easier to calculate the power spectrum than other 
resummation methods even in the presence of halo bias and 
redshift-space distortions. The computational cost is similar 
to that of the SPT. 

We focus not on power spectrum but on two-point correla- 
tion function, because we do not suffer from shot noise effect 
in the correlation function. While 2-loop corrections in the 
iPT generally might have an impact on the correlation func- 
tion on quasi-linear scales [11411 . we use 1-loop iPT for sim- 
plicity in this paper 

This paper is organized as foUows. We first review the theo- 
retical predictions of power spectrum and correlation function 
in Section [III We describe the details of A-body simulations 
and a method to compute the correlation functions of mat- 
ter and halos from A-body simulations in Section |III1 After 



showing the resuhs of matter and halo correlation functions 
and its cross-correlation function in Section |IV] we then show 
the main results of this paper in Section|V] Finally, SectionlVTl 
is devoted to our conclusion. 



II. ANALYTIC MODELS 

In this section, we briefly review two theoretical models: 
the iPT with nonlocal bias model and the SPT with local bias 
model, which are compared with A^-body simulation results. 



A. Predictions of integrated Perturbation Tlieory 



We use iPT lllOll to investigate how the cross-correlation 
coefficient behaves on quasi-linear scales for various halo 
masses and redshifts. It is convenient to write down the power 
spectrum predictions of the iPT based on multipoint propa- 
gators recently introduced in Bemardeau et al. llisll . Using 
the multipoint propagators r^"\ one-loop power spectrum be- 
tween object X and Y based on the iPT can be written as (the 
full derivation is given in uM) 
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where indices X and Y are either matter 'm' or halo 'h' in this paper, Pi,{k) is the linear matter power spectrum, and vertex factor 
n{k) is given by 



n{k) = exp 
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The normalized multipoint propagators T-'^ and T^^' are given by 
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where indices 'm' and 'h' denote the matter and halo, 
and we assume the second-order renormalized bias function 
c|j (fci, fca) depends only on the magnitudes of wavevectors, 
k\ = \ki\ and k2 = Ifcal- From Equations (|3ll-(|6|l, we can 
easily understand that the matter result is recovered when 
ci = Cu - Q- Here c. and c. are renormalized bias func- 

A li 

tions in Lagrangian space introduced in Matsubara IllOll and 
obtained as 
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for a mass range Mi < M < M2 (see Equations 64 and 108 
of Matsubara uM ). where v = 6c/a-(M) is a function of mass 
M, and 6c is the critical overdensity for spherical collapse. 
In an Einstein-de Sitter cosmology, the critical overdensity 
is dc ~ 1 .686, while it shows weak dependence on cosmol- 
ogy and redshifts in general cosmology lll8lll9ll . and thus we 
use Henry 11911 fitting formula to include cosmological depen- 



dence. The function cr(M) is the root-mean-square linear den- 
sity field smoothed with a top hat filter of radius R enclosing 
an average mass M = poA-nR^ /3; 
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where po is the mean matter density of the universe and ji (x) 
is the first-order spherical Bessel function, /mf(t^) is the scaled 
differential mass function defined as 1I20I1 
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where n(M) is the comoving number density of halos with 
mass M. The quantity /mf(v) is frequently used in the lit- 
erature and there have been several analytic predictions lEll - 
I23II and fitting formulas [e.g.. l20ll24l - l29ll . In this paper, we 



use Bhattacharya et al. II29I1 fitting formula for mass function 
which shows better agreement with our simulations ||30|. c^ 
is given as (see. Equations 92, 95, and 96 in Matsubara lUTIl ') 
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with 



Ao(M) = 1, 

Ai(M) = 1 +6cb\(M), 



(12) 
(13) 



where b]^ is the Lagrangian bias function for the halo bias. 

Theoretical two-point correlation function can be expressed 
in terms of the power spectrum as 
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where F2 is the 2nd-order mode coupling kernel in SPT: 

5 Iki-kolki ko\ 2(ki-k2\^ 

Taking Fourier transforms, we then obtain corresponding 
correlation functions given by 
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^hm('-) = ^i^nl('-) + b2A(r), 
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where ^nl is the nonlinear matter correlation function, A(r) 
and B(r) are the Fourier transform of A(k) and B(k). Note 
that B{r) - ^\{r) - crl Soir) where ^l('") is the linear matter 
correlation function, crl - J d^q P^(q)/(2n)^, and doir) is the 
Dirac delta function. 



III. W-BODY SIMULATIONS 



B. Standard perturbation theory with local bias model 

In the SPT formalism, we consider the local determinis- 
tic nonlinear biasing model. Following Fry and Gaztanaga 
Oin . we restrict the consideration on large scales in Eulerian 
space and assume that the halo density can be described by 
a smoothed function !F((5m) that depends only on the matter 
density. We can expand 'F in a Taylor series around 6m such 
that 
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where (5^ is the nonlinear matter density. We then combine 
this expansion with SPT, which expands the matter density 
perturbations into a series 6^ - ^m + ^L + ■ ■ • > where SiJ 
is the linear density field and 6^ is of order [^Jn ]". At the 
next-to-leading order, we can obtain the auto- and cross-power 
spectrum of halos as lllSl [3211 

Phh(^) = biPNL(k) + 2bib2A(k) + -j-B{k) + N, (16) 
P^^mik) ^ byPni^ik) + b2A{k), (17) 

where b\ and b2 are the renormalized bias parameters, A^ is the 
renormalized shot noise, and Pnl(^) is the nonlinear matter 
power spectrum. b\ and b2 should be determined empirically 
or treated as free parameters. In this paper, we will exam- 
ine both cases in Section |Vl The terms A{k) and B{k) can be 
obtained as 
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A. Simulation parameters 

To obtain accurate predictions of the cross-correlation co- 
efficient, we resort to use the high-resolution A^-body simu- 
lations of structure formation. To perform the A^-body simu- 
lations, we use a publicly available tree-particle mesh code. 
Gadget! iH]. We adopt the standard ACDM model with the 
matter density Qm = 0.265, the baryon density Qb — 0.0448, 
the dark energy density Qa = 0.735 with equation of state pa- 
rameter w = -1, the spectral index n, = 0.963, the variance 
of the density perturbations in a sphere of radius 8/j"'Mpc 
0-% - 0.80, and the Hubble parameter /; = 0.71. These cos- 
mological parameters are consistent with the Wilkinson Mi- 
crowave Anisotropy Probe 7yr results [J34J. We performed 
two types of simulations both with Np - 1024'' particles in 
cubic boxes. First type has a side Lbox = lOOO/j^'Mpc with 
softening length r, being 50/j"'kpc, and the second type has 
a side Lbox = 2000/i"'Mpc with softening length r, being 
lOO/j^'kpc. These two types are named as LIOOO and L2000, 
respectively. The initial conditions are generated based on the 
2nd-order Lagrangian perturbation theory [351 36| with the 
initial linear power spectrum calculated by CAMB [37]. The 
initial redshift is set to Zim - 36 for LIOOO and Zim =31 for 
L2000. We perform A^run = 30 and 10 reahzations for LIOOO 
and L2000, respectively. We summarize simulation parame- 
ters in Table U The LIOOO simulations used in this paper are 
the same as LIOOO used in Sato and Matsubara IBOll . 

We store outputs at z = 3.0, 2.0, 1.0, 0.5, 0.3, and 
and identify halos for each output using a Friends-of -Friends 
(FOF) group finder with linking length of 0.2 times the mean 
separation IBSIl . We select halos in which the number of par- 
ticles, Np, is equal to or larger than 20 which corresponds 
to the halos with masses 1.37 x IO'^/i^'Mq for LIOOO and 
1.10 X IO'^^/i^'Mq for L2000. Then we divide halos into five 
mass bins to keep track of their different clustering properties. 
The average number and mass of halos among realizations for 
redshifts are listed in Table |ll] The halo catalogs of Bin 4 in 



TABLE I: Parameters in high- and low-resolution A'-body simulations: the matter density Q^,, the dark energy density Qa, the baryon density 
Q.\,, the Hubble parameter h, the spectral index ;ij, the variance of the density perturbations at 8/j"'Mpc erg, the box size L^ox, the number of 
particles A',,, the initial redshift Zm\, the softening length r„ and the number of realizations N^„. 



Name 



iZnn 



Ha 



Oh 



OS 



N„ 



A'm 



LI 000 (high resolution) 
L2000 (low resolution) 



0.265 0.735 0.0448 0.71 0.963 0.80 lOOO/i'Mpc 1024^ 36 50/r'kpc 30 
0.265 0.735 0.0448 0.71 0.963 0.80 2000r'Mpc 1024^ 31 lOOr'kpc 10 



TABLE II: Properties of halo catalogs of high- and low-resolution A'-body simulations for each mass bin. We use the 5 halo catalogs 
abbreviated to "Bin 1",. . . , "Bin 5". ^h and Mj, are the average halo numbers and average halo masses at various redshifts. 
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L2000 
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L2000 is constructed so that the halo mass range is the same 
as that of Bin 3 in LIOOO, as shown in Table HIl Since volume 
of L2000 simulations is bigger than that of LIOOO simulations, 
number of halos with a certain mass are larger for L2000 sim- 
ulations. 



B. Analysis: Two-point correlation functions 

To calculate the two-point correlation function of dark mat- 
ters from A'-body simulations, we first randomly choose num- 
ber of particles A';,,,- = 196^ and 256^ for LIOOO and L2000. 
For dark matter halos, we use all halos in each bin. Then we 
directly count the A'-body particle and/or halos to calculate 
two-point correlation function instead of using method with 
Fast Fourier Transform llsoll . We choose r, to be the center of 



The shot noise corrections in the halo power spectrum are 
subtle. If the dark matter halos are regarded as a Poisson pro- 
cess, we can easily subtract shot noise effect using number 
density of halos «/,. However Smith et al. [i39i] found that 
this standard correction method is not exactly correct for ha- 
los, particularly for those of large mass. This is probably be- 
cause in order to identify halos using the FOF algorithm, we 
automatically impose that distances between halos are larger 
than the sum of their radii, or they would have been linked as 
bigger halos. Thus, the shot noise effect is scale dependent 
and it is difficult to correctly subtract the effect of shot noise. 
Therefore, we use the coiTelation function instead of using the 
power spectrum, because the shot noise effect in the correla- 
tion function is weaker than that in the power spectrum. 



the /-th bin, i.e., n = (r™'" + r"")/2, where r™" and r™'"' are 
the minimum and maximum distances of the /-th bin. 



z=0 
M^ = (1.37-4.11)xlO'2[h-'Mg] 




cross-correlation coefficient between matter and halo density 
fields, defined as 



10 so 30 40 50 60 70 80 90 100 110 120 
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FIG. 1: Top panel: Congelation functions of matter, halo, and their 
cross-coiTelation function at redshift z = 0, multiplied by a factor of 
1^. For halo mass range, we consider 1.37 < Mh/(I0'^/7"'Mo) < 4.11 
(Bin 1). The symbols are the results obtained from A'-body simula- 
tions. The solid lines are the results of integrated perturbation theoiy 
(iPT) llOll . Bottom panel: Ratio of correlation functions measured 
from A'-body simulations to those from the iPT. 



IV. CORRELATION FUNCTIONS 

Before presenting the results of cross-correlation coeffi- 
cient, we compare the A'-body simulation results with the iPT 
for correlation functions themselves. 

Figure [1] shows the results for the correlation functions of 
matter and halos, and their cross-correlation function at z = 0. 
We use halo catalog of "Bin 1" shown in Table U The am- 
plitude of halo-halo correlation is smaller than that of matter- 
matter correlation, because the halo bias b^ in this halo range 
is 0.904 (less than 1). The error bars describe the 1-cr error 
on the mean values obtained from 30 realizations. The error 
bars increase on large scales because of the finite size of the 
simulation box. The iPT predictions agree with A'-body simu- 
lation results down to r ~ 25/?"' Mpc within a few percent for 
all correlations. In Section |V] we will see that a range of a 
few percent-level agreement in cross-correlation coefficient is 
extended more than that in the correlation functions. 



V. CROSS-CORRELATION COEFFICIENT 

In the framework of local biasing model, density field of 
galaxies and their halos should be a stochastic function of the 
underlying dark matter density field [J9J]. The stochasticity is 
very weak on large scales, while it becomes more important 
on small scales ||40| - |42I1 . 

One of the characteristic parameters of stochasticity is the 
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where ^mm('"), ^hh('"), and ^hm('") are the matter and halo auto- 
correlation functions, and their cross-correlation function, re- 
spectively. The cross-correlation is the measure of the sta- 
tistical coherence of the two fields ||431 - |47I1 . If any scale- 
dependent, deterministic, linear-bias model is assumed, we 
have Tec = 1 . Therefore, deviations of cross-correlation co- 
efficient from unity would arise due to both the nonlinearity 
and stochasticity of bias. 

Figure [2] shows the cross-correlation coefficient between 
matter and halo density fields at z = 0, 0.3, 0.5, 1.0, 
2.0 and 3.0. The cross, square, and triangle symbols are 
ffie A'-body simulation results measured from 30 realizations 
for halo masses 1.37 < Mh/(10'VMa) < 4.11 (Bin 1), 
4.11 < Mh/(10'-/i-'Mo) < 12.32 (Bin 2), and 1.23 < 
Mh/(10"h'^Mo) < 3.70 (Bin 3), respectively The error bars 
describe the l-o" error on the mean value obtained from 30 
realizations. We do not plot the results in which sum of the 
1-cr error bars in a range of 5 < r,/(/!"'Mpc) < 100 is larger 
than 0.12, i.e., 25<r,<iooO"' > 0.12. It should be noted that 
halos in each bin are more biased as increasing redshift, be- 
cause we impose the same halo mass ranges for each bin. The 
solid curves show the iPT predictions. The iPT obtains good 
agreements with simulation results down to r ~ 15/2" 'Mpc 
within a range of error bars for all redshifts and halo mass 
ranges we have considered. Particularly at z = 3.0, the iPT 
well reproduces the simulation result down to r ~ 6/z 'Mpc. 
Below, we will see fractional differences between simulation 
results and theoretical predictions in Figure |5] to discuss the 
percentage error The difference between the iPT and simula- 
tion results on small scales probably comes from the fact that 
the iPT breaks down on small scales (see. Figure [T} OOi 14811 . 
One can see that the iPT prediction on small scales are almost 
flat unlike the simulation results. This is probably because the 
asymptotic behaviors of the correlation functions based on the 
iPT are almost the same each other (see, Figure[T}, and the iPT 
should not be applied on such small scales at any rate. 

We also plot a simple model derived from Equations (ISTt 
and (|22] | as dotted curves and it is expressed as ll32ll 
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by using following approximations; A(r) « ^Nhir) and 
B{r) « ^NL(r). 

To empirically estimate bi and b2, we use general relations 
between local bias parameters in Lagrangian space and Eu- 
lerian space, which are derived in spherical collapse model 
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FIG. 2: The cross-correlation coefficient between matter and halo density fields at z = 0, 0.3, 0.5, 1.0, 2.0, and 3.0. For halo mass ranges, we 
consider 1.37 < Mh/(10''/i-'Mo) < 4.11 (Bin 1), 4.11 < Mh/ClO'^/^-'Mo) < 12.32 (Bin 2), and 1.23 < Mh/(10"/i-'Mo) < 3.70 (Bin 3). We 
do not plot the results in which sum of the 1-cr error bars in a range of 5 < r,7(/j"'Mpc) < 100 is larger than 0.12, i.e., 25<;<ioo^i > 0.12. The 
symbols are the results measured from A'-body simulations. The solid lines are the results of integrated perturbation theory (iPT) (kJ while 
the dotted lines are the results of simple model derived from standard perturbation theory with local bias model fEquation l24t . To empirically 
estimate bi and ^2, we use the relation in Equations ( 125 l l and ( I26l l and then simply substitute b^ and fc| with hi and b2. 
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FIG. 3: Same as Figure|2] but for the dotted lines, we fit b2/bi to the A'-body simulation results using chi-square fit. 
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FIG. 4: The cross-correlation coefficients with halo mass ranges: 1.23 < Mh/(10"/;"'Mo) < 3.70 (Bin 4), and 3.70 < Mh/(10"/r'Mo) < 
11.09 (Bin 5), given at redshift z = 0, 0.3, 0.5, and 1.0. We do not show the results in which sum of the l-cr error bars in a range of 
5 < rj/(h^'Mpc) < 100 is larger than 0. 12. The tr iang le and circle symbols are the results of A'-body simulations. The solid lines correspond 
to results of integrated perturbation theory (iPT) lUoll while the dotted lines correspond to results of standard perturbation theory with fitted 
bias model (Equation l24b . 



where b^ and i>| are Eulerian bias parameters. Note that both 
the Eulerian bias parameters b^ and the Lagrangian bias pa- 
rameters b}; are local and independent on scales. In this phe- 
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and b2- To calculate b\ and bl;, we use II12I1 
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for halos in a mass range M\ < M < Mo- The simple model 
(Equation [24l i with the above estimates of bias parameters 
shows better agreement with simulations for higher redshifts 
(i.e., more biased halos). We find that the cross-correlation 
coefficients of halos with b\ > 2 are well described in this 
method over all scales we considered. For lower redshifts 
(i.e., less biased halos), the simple model deviates more from 
the simulation results. 
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FIG. 5: Fractional differences between A'-body results and theoretical predictions are shown in percents, at redshifts z = 0, 0.3, 0.5, 1.0, 2.0, 
and 3.0, from bottom to top. Left panels: iPT, middle panels: SPT, and right panels: SPT with fitting. The red, blue, and green symbols are 
the results of Bin 1, 2, and 3. 



symbols are the same as the Figure |2] Fittings are done in a 
range of 5 < r/(/!"'Mpc) < 70. In the fitting case, an im- 
provement from the above empirical method is little for cases 
of high bias, but is important for cases of low bias. The simple 
model with fitted bias replicates the simulation results over all 
scales at 1 < z < 3. We can see that the cross-correlation 
coefficients estimated from A^-body simulations have compli- 
cated behaviors in quasi-linear regimes at low redshifts, which 
cannot be described in the simple model. We will describe 
percentage error in Figure |5]later. 



Figure |4] shows the results for the cross-correlation 
coefficient of large halos with mass ranges: 
1.23 < Mh/(10'^r'Ma) < 3.70 (Bin 4) and 
3.70 < Mh/ClO'^/i-'Mo) < 11.09 (Bin 5) at redshift 
z = 0, 0.3, 0.5, and 1.0. The triangle and circle symbols 
are the simulation results of Bin 4 and Bin 5 estimated from 
10 realizations of L2000. The soUd and dotted lines are the 
predictions of the iPT and SPT with fitted bias, respectively. 
As in Figures |2] and |3] the iPT shows nice agreement with the 
simulation results on large scales even in large halo masses. 



10 



&^ 



- 



-1 



-1 - 



-1 - 



-2 



-r-r-q 

z=1.0 



-1 1 1 — I — r- 



Jl 



U .. , ti 



,/''«I-'f 



li 



-^ll^ 



■li' 



I I I I 



z = 0.5 



H h 



¥\i 



H 






<-+-H 1 \ — I I J 1 



z = 0.3 



^^Aifif 



^ 



z = 



H \ — I I J 






iPT 



-i-n-f 

z=1.0 



I.I 



..:.::.yftwiii 



^h' 



1 1 1 1 



H \ — I— H 



' = °-' * ] A ill] I 

■■■■mJm 



I I I I I \ \ — M l!! 

z = 0.3 



#^ 



ih. 



_L 



_L 









,'iL 



l!i! 



^ 



H- 



H h 



z = 









SPT 




J I I I I I 



--1-0 llli IT 



tjg.^j 



-I-+-H 



H ! — Hr-f 






! I I I 



H ! — I l -J 1. 



,- z=i 






ill 



H ! — I MtI 



z=0 



"t Tifi^WlIlIlT!! 



_L 



SPT w fit 

J I L_L_ 



5 10 50 100 10 50 100 10 50 100 

r[h"' Mpc] r[h"' Mpc] r[h"' Mpc] 



FIG. 6; Same as Figure|5] but the results of Bin 4 (green triangle) and 5 (magenta circle) at z = 0, 0.3, 0.5 and 1.0. 



The simple model with fitting also reproduces the simulation 
results for large halo masses. However, fitting values of bi/bx 
are not, in general, the same as those obtained from other 
statistics such as power spectrum and bispectrum, because bi 
and i>2 are renormalized. 

To clearify how well theoretical models predict the A^- 
body results, we plot fractional differences between A^-body 
simulation results and theoretical predictions, [r^csimif) " 
rcc.theoif)]/ fcc.theoif), as shown in Figures |5] and |6] These Fig- 
ures show that the iPT agrees with simulation results down to 
r ~ 15 (10)/!"'Mpc within 0.5 (1.0) % for all redshifts and 
halo masses we considered. It should be noted that the iPT 



does not have any fitting parameter The SPT with empiri- 
cally determined bias reproduces A^-body simulation results 
down to r ~ IO/2 'Mpc within a percent-level for all redshifts 
except for z = (see. Figure |5]l. In the SPT with bias deter- 
mined by fitting, a percent-level agreement is achieved over 
wide separation angles for all redshifts. However, the fitted 
parameters bi and b2 are different from b'^ and b^ which can 
be determined by other method, e.g., power spectrum and bis- 
pectrum. 
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VI. CONCLUSION 

In this paper, we have used 40 large cosmological A^- 
body simulations of the standard ACDM cosmology to in- 
vestigate the cross-correlation coefficient between halo and 
matter density fields over a wide redshift range. The cross- 
correlation coefficient is crucial to extract information of 
the matter density field by combining galaxy clustering and 
galaxy-galaxy-lensing measurements. Since the first attempt 
to detect galaxy-galaxy lensing 14911 . its ability to constrain 
cosmological parameters has been shown ISft. 

We compare the simulation results with theoretical predic- 
tions of the iPT and simle models of bias with SPT. The iPT 
predicts the simulation results down to r ~ 15 (10)/i"'Mpc 
within 0.5 (1.0) % for all redshifts and halo masses we con- 
sidered. To improve the prediction, the 2-loop correction to 
the iPT might be important. In the SPT with local bias model, 
bias parameters are renormalized and therefore they are de- 
termined empirically or treated as free parameters. The SPT 
with empirically determined biases with spherical collapse 
model shows better agreement with simulations for more bi- 
ased halos on small scales, although this model does not re- 
produce complicated behaviors of the simulation results on 
quasi-linear scales at low redshifts. The SPT with biases 
determined by fitting improves the predictions but the situ- 
ation is almost the same at low redshift. Thus, the iPT ac- 
curately predicts the cross-correlation coefficient as long as 
quasi-linear scales are considered. 

Let us finally comment on convolution Lagrangian pertur- 
bation theory (CLPT), which is recently proposed by Carlson 



et al. 114811 . The CLPT applies additional resummations on top 
of the simple LRT (restricted iPT with local Lagrangian bias), 
and its prediction significantly improves the simple LRT for 
correlation function in real and redshift spaces on small scales. 
Therefore, it might be possible that the CLPT gives a better 
prediction for the cross-correlation coefficient between mass 
and halos and agrees with simulation results on small scales. 
Although it is important to examine how well the CLPT pre- 
dicts, we leave it for future work. 

In this paper we focus on fundamental features of bias 
stochasticity by the methods of numerical simulations and the- 
oretical models. We believe the results of this paper could be 
a crucial first step to understand the galaxy biasing for future 
precision cosmology. 
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